Phase transitions and perfectness of fluids in weakly coupled real scalar field theories 
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We calculate the ratio rj/s, the shear viscosity (n) to entropy density (s), which characterizes 
how perfect a fluid is, in weakly coupled real scalar field theories with different types of phase 
transitions. The mean-field results of the n/s behaviors agree with the empirical observations in 
atomic and molecular systems such as H2O, He, N, and all the matters with data available in the 
NIST database. These behaviors are expected to be the same in N component scalar theories with 
an O(N) symmetry. We speculate these n/s behaviors are general properties of fluid shared by 
QCD and cold atoms. Finally, we clarify some issues regarding counterexamples of the conjectured 
universal bound n/s > l/4ir found in Refs. [l6l.[l7l]. 
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Quantum chromodynamics (QCD) is believed to un- 
dergo rapid transitions from hadronic phases to a quark- 
gluon plasma (QGP) phase at high temperature T and 
to quark matter phases at high quark chemical potential 
ji (see P, 0, E| for reviews). Lattice results show that 
the phase transition of the hadronic matter to QGP at 
finite T with \x = is likely a crossover At finite jx 
and T = 0, there is no reliable lattice result due to the 
severe fermion sign problem. However, arguments based 
on a variety of models show that the phase transition is 
of first-order. This first-order phase transition turns into 
a crossover at smaller ji and finite T at the QCD critical 
end-point (CEP) @. There are lots of interests in the 
CEP. Recently, it was proposed to probe the CEP by us- 
ing the ratio of shear viscosity 77 to the entropy density s 
of QCD @. 

Shear viscosity r\ characterizes how strongly particles 
interact and move collectively in a many-body system. 
In general, the stronger the interparticle interaction, the 
smaller the shear viscosity (here r\ is normalized by the 
density). It was conjectured Q that no matter how 
strong the particle interaction is, rj/s has a universal min- 
imum bound 1/47T in any system. This bound is moti- 
vated by the uncertainty principle and is found to be sat- 
urated for a large class of strongly interacting quantum 
field theories whose dual descriptions in string theory in- 
volve black holes in anti-de Sitter space 0j- 

There are two important questions regarding this con- 
jecture. The first one is: Is this rj/s bound truly uni- 
versal? By definition, there is no proof of this conjec- 
ture yet. From experimental observations, the bound is 
well satisfied in matters like H2O, N and superfiuid He 
Q. For cold fermionic atoms with an infinite scatter- 
ing length (the unitarity limit), the bound is satisfied 
but rj/s is close to the bound near the phase transi- 
tion temperature T r . [fill . Similar behavior is found for 
QCD at u, = M. | l ( J . Relativistic heavy ion collisions 
(RHIC) [Tl|, [STlirill and lattice computations of a 
gluon plasma suggest that rj/s of QCD is close to 
the minimum bound at just above T c . But one cannot 
conclude whether the bound is violated based on current 



precision. On the other hand, interesting counterexam- 
ples of the bound have been constructed 
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[T7I ] and a 

possible modification of the bound due to some string 
excitations in the dual theory [ID, 0, [2(| • 

The second question is: What is the general behavior 
of rj/s, especially for QCD and cold atom systems which 
arc of high interests but not well known? The answer to 
this question is interesting on its own. It also has inter- 
esting applications. For example, locating QCD CEP by 
•q/s requires this information. Also, knowing when a sys- 
tem will reach its minimum rj/s is important for testing 
the minimum bound conjecture. By now it is known that 
several different systems have qualitatively the same rj/s 
behaviors. For example, H2O, N, and He all have the 
minimum rj/s near the liquid-gas phase transition tem- 
perature. The same is true with QCD at fj, = J^, [l(| 
and near the nuclear liquid-gas phase transition [2l|. It 
is also true in cold atoms in the unitarity limit [8|. It 
would certainly be interesting to further map out the de- 
tailed rj/s structure in the phase diagram of QCD and 
other systems. 

In this letter, we address the above two questions at 
the same time. We study how rj/s behaves in the simplest 
field theory — a real scalar field theory. The resulting rj/s 
in first-, second-order phase transitions and crossover be- 
haves the same way as in H2O, N, He, and all the matters 
with data available in the NIST database @, HH, Hi] . We 
argue that this agreement might hold when the theory is 
generalized to N components with an O(N) symmetry, 
which is the low energy effective field theory of a big class 
of systems. We then speculate these behaviors might be 
general properties of fluid which are shared by QCD and 
cold atoms. Furthermore, this simple theory allows us to 
clarify some issues regarding the counterexamples of the 
rj/s bound [l6L in a more transparent way. 

We will study a real scalar theory with the Lagrangian 

C=\{d,cpf- l -a ( j ) 2 -\b^- l -c^. (1) 

This theory is invariant under (f> — > — (f> and has a Z2 
symmetry. There could be two additional terms with 
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dimension six: 3 <9 2 0, <j>d 2 d 2 4> (the other terms are re- 
lated to these ones by integration by parts). These terms 
can be removed by field redefinition or, equivalcntly, by 
applying the equation of motion. The inclusion of the 
dimension six terms shows that this is an effective field 
theory, which is valid under the cut-off scale and is 

renormalized order by order in the momentum expansion 
Py/c, p being a typical momentum scale in the problem, 
a, 6, and c are renormalized quantities and the countert- 
erm Lagrangian is not shown. The renormalization con- 
dition is that the counterterms do not change the particle 
mass and the four- and six-point couplings at threshold. 
We will discuss the following cases: l)c = 0,£>>0, 
a > 0, the system is always in the symmetric phase. 2) 
c = 0, 6 > 0, a < 0, the vacuum at T = breaks the Zi 
symmetry spontaneously. However, the symmetry is re- 
stored at higher T with a second-order phase transition. 
3) Adding an explicit symmetry breaking term 8L = H(f) 
to the Lagrangian of 2) to model a crossover. 4) c > 0, 
6 < 0, a > 0, the broken symmetry is restored at high T 
with a first-order phase transition. 

We will focus on the case of weak coupling and com- 
pute the mean-field effective potential via the standard 
Cornwall-Jackiw-Tomboulis (CJT) formalism [23[ which 
has the one-particle irreducible diagrams included self- 
consistcntly. The effective potential in the CJT formal- 
ism reads 
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where U{<j>) — a/2 cj) 2 + 6/4 4 + c/6 4> & is the tree- level 
potential, and S(So) is the full(tree-level) propagator: 



S-\K, 



-K 2 + m 2 $) ,SoHK, 



with the tree-level mass m 2 . 



-K 2 + m 2 S) , 
(3) 
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In this work, we neglect the non-tadpole type loop 
diagrams. This "Hartree approximation" is good when 
T > \b\ 1/2 T c for c = and when T > |a| 1/2 for c > 0. 
Thus, to the order we are working, the 2 PI potential Vi 
only includes 
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where L (0) = J K S(K,<j>). The self-consistent one- and 
two-point Green's functions satisfy 
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<i>=4>o,s=S(4>o) 
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' ss 







(5) 



<t>=<t>o,S=S(tfro) 



This allows us to solve (j} Q and m through the coupled 
equations: 

<t> (a + beg + c<t>t + (36 + lOc0g)£(0o) + 15cL(c/> ) 2 ) = . 
m 2 ~m 2 = 3(6 + lOc^)i(0 o ) + 15cL(0 o ) 2 . (6) 
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FIG. 1: rj/s vs. T for cases with a second-order phase transi- 
tion (solid curve), a crossover (dash-dotted curve), and with 
no phase transition for massive field (dashed curve) and mass- 
less field (dotted curve). Parameters can be in arbitrary units. 



The entropy density of the system is given as s = 
—dV(4>o)/dT, while the shear viscosity r\ is calculated us- 
ing the Boltzmann equation. It is proven that in a weakly 
coupled scalar field theory with quartic and cubic terms, 
summing the leading order diagrams for 77 is equivalent 
to solving the Boltzmann equation with effective T de- 
pendent mass and scattering amplitudes [25|. Thus, one 
can directly apply the Boltzmann equation to compute 77 
for cases with c = for both the symmetric (a > 0) and 
symmetry breakin g (a < 0) cases. Furthermore, since 
the proof of Ref. [25| does not use properties that are 
restricted to scalar theories, the conclusion is expected 
to hold for more general theories with weak couplings, 
including QCD in the perturbative regime (26[. Here, we 
also apply it to the c > case with a first-order phase 
transition. 

The two-particle clastic scattering amplitude, which 
governs particle collisions in the Boltzmann equation, is 



iT = A 4 + A 2 
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where s, t and u are Mandelstam variables, and A 3 = 
6<M6 + + 1OcL(0o)) and A 4 = 12(§ + 5c0g + 

hcL (<f>a)) are effective couplings. 

Now we discuss the behavior of 77/s, starting from the 
simplest case: a = c = and b > 0. In this theory, 
the only scale in the problem is T and the system is 
always in the symmetric phase. Thus, both r\ and s are 
proportional to T 3 on dimensional ground and ij/s is T 
independent (up to the logarithmic running of b which 
we neglected). 

When the mass term is added, 77/s is no longer a con- 
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FIG. 2: rj/s vs. T for the ground state (solid curve) and 
metastable state (dash-dotted curve) in a first-order phase 
transition. Parameters can be in arbitrary units. 

stant. rj/s curves for b > 0, c = and (i) a > (ii) 
a < arc shown in Fig. 1. (Wc have used b = 0.1. The 
expansion parameter is ~ A/(4tt) 2 , where A = 66. (27j ) 
For (i), there is no phase transition, rj/s is always mono- 
tonically decreasing. For (ii), rj/s reaches its minimum 
and develops a cusp at the T c of the second-order phase 
transition. There is clearly a qualitative difference in the 
rj/s behavior between cases with and without a phase 
transition. In fact, in the high T expansion (bT 2 ^> \a\) , 

with k = 192 + 4S)Vb + 0(b), consistent with rj com- 
puted in [25|,[28|. Thus, the sign of a determines whether 
rj/s is decreasing or increasing at high T. As T — ► 0, s 
approaches zero exponentially (the excitations are mas- 
sive) while r\ approaches zero by power laws. Thus, rj/s 
is decreasing in both cases at low T. (This feature is 
not affected by the Hartree approximation used at low 
T.) One concludes that rj/s is not monotonic in T for 
a second-order phase transition. Although our mean- 
field result would be modified very close to T c (with 
\T — T c \ /T c < 0(b)), we argue that the most natural 
scenario for rj/s in a second-order phase transition is to 
have a single local minimum at T c (where m also reaches 
its minimum, m = 0), because there is no other T that 
is more special than the others to develop another mini- 
mum in rj/s. Analogously, the most natural scenario for 
rj/s with no phase transition is monotonic decreasing be- 
cause no T is more special than the others to develop a 
minimum. Finally, for a second-order phase transition, 
one naturally expects a cusp at T c , while the cusp is 
smoothed out in a crossover (H ^ 0). These "natural- 



ness" arguments do not depend on whether the coupling 
is weak. In fact, there is no b that is more special than 
the others such that one does not expect the rj/s behavior 
to change in the strong coupling cases (|6| ^> 1) either. 

For the first-order phase transition (a > 0, b < 0, and 
c > 0), the result for rj/s is shown in Fig. 2. There 
is a discontinuity at T c , as expected. The rj/s is con- 
cave in general and its minimum is reached at T c from 
the low T side. The qualitative behaviors above and be- 
low T c are similar to those of a second-order phase tran- 
sition. The main difference is the discontinuity. It is 
not intuitively obvious why this should be the case be- 
cause even though the effective potentials of the first- and 
second-order phase transitions look similar near the min- 
ima (thus the masses behave similarly) the Lagrangians 
are still quite different. It turns out this is simply be- 
cause the scattering amplitude of Eq. (|7|) is dominated by 
the A4 term in both cases as long as T 2 > \a\. Also, 
the fact that the minimum of rj/s is reached at T c from 
below, instead of above, is physical but nontrivial. It is 
physical, because particle interaction is stronger in the 
low temperature side. It is nontrivial, because the dis- 
continuities in rj and s are of the same sign. So it is not 
obvious the discontinuity of rj/s has to be governed by rj. 

The above-mentioned behaviors of rj/s in first-, second- 
order phase transitions and crossover are the same as 
those in H2O, N, He, and all the matters with data avail- 
able in the NIST database @, SH, |||. We expect the 
agreement still holds when the scalar field is generalized 
to N components and the Lagrangian has an O(N) sym- 
metry. This is because the arguments used to explain the 
rj/s behaviors do not depend on N, except for Eq.® in 
the c = case. However, in the O(N) theory, the pref- 
actors in Eq.© change [2!| but the signs remain. Thus, 
the arguments still hold. 

Given the general agreement in the rj/s behaviors of 
H2O, He, N, O(N) scalar field theories, and QCD and 
cold atoms in some known cases, it is conceivable that the 
rj/s behaviors shown in Figs. 1 and 2 are general prop- 
erties of a large class of fluids, including QCD and cold 
atoms. If this is correct, then the best place to measure 
the minimum rj/s to test the minimum bound conjecture 
is near T c . In particular, for a pure gluon theory which 
has a first-order phase transition, the minimum rj/s could 
be just below T c . This strongly motivates future lattice 
computations to be carried out in the lower T regime. In 
addition, this also provides theoretical support for locat- 
ing QCD CEP by measuring rj/s. There is some subtlety 
regarding the CEP. In H2O, He, or N, a cusp in rj/s is 
observed at the CEP with experimental resolution. The- 
oretically, it is argued that 77 has a mild divergence (and 
hence rj/s) at QCD CEP [301 as in the Model H accord- 
ing to Hohenberg and Halperin's classification [3l| (but 
see [HI for a caveat). Experimentally, divergence of rj in 
3 He liquid-gas transition critical point was observed [33j 
but not for the tricritical point of 3 He- 4 He binary fluid 
[34| . Both cases belong to Model H and we re p redicted 
to have the same weak divergence by Ref. [3l[. In any 
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case, it is unlikely that RHIC has the sensitivity to this 
effect. 

Finally, we comment on the counterexamples of the 
minimum 77/s bound conjecture constructed in (l6l. Il7l]. 
A common feature of those counterexamples is a large 
number of flavor g. In the weak coupling and small den- 
sity limits, one can compute 77 (also using the Boltzmann 
equation) and s reliably. In the non-relativistic exam- 
ples employed in pH , [l7| , particle numbers are fixed or 
controlled by chemical potentials. It was found that 77 
does not scale with g (because it is normalized to the 
density already) but s does. Thus, 77/s oc 1/ {Log(g)a 2 a ) 
as g — > 00 (a s is the two-body scattering length) and the 
1/4-7T bound was violated. There is an important differ- 
ence in those counterexamples on whether flavor chang- 
ing is allowed during the collisions. In the O(N) scalar 
field theory <fii<fii — ► 4>i<j)i (i = 1,2...N) collision can 
happen. Thus, it is a model with flavor changing. To 
the leading order in the b expansion, 77/s oc l/(Nb 2 ) for 
a = c = 0, similar to Eq.© but with an extra factor of N 
from s. Naively the conjecture is violated when N — > 00. 
However, the expansion parameter is (Nb). One needs 
b ex 1/N or smaller (and hence 77/s cx N or larger) to be 
consistent with the weak coupling treatment. Thus, one 
cannot conclude that the O(N) scalar field theory vio- 
lates the bound in the large iV limit based on weak cou- 
pling argument. The same is true to the non-relativistic 
flavoring changing counterexamples of [l6l Il7| (a s oc 1 / g 
is required). In the strong coupling regime, the b oc 1/N 
scaling seems not necessary. However, without the scal- 
ing, the P function of b blows up in the large N limit. It is 
not clear whether a sensible non-pcrturbativc theory can 
be constructed this way. Our example sets no constraint 
on flavor conserving cases in [h| [lj], though. 

In Ref. [TU, there is discussion about whether the 
77/s bound conjecture should be applied to metastable 



states. In the definition using Kubo's formula, 77 is re- 
lated to the linear response of an ensemble. It is hard to 
exclude the ground state from the ensemble to compute 
the metastable state 77. A similar problem happens in 
computing 77 with the Boltzmann equation in which the 
dissipation of perturbations away from thermal equilib- 
rium gives rise to 77. The ground state information is en- 
coded in the equilibrium distribution. However, experi- 
mentally it makes perfect sense to measure 77 of a material 
as long as it is stable during the time of the measurement. 
Water is a sharp example of this [l|| . In electroweak in- 
teraction, water is metastable to proton fusion of the two 
hydrogen atoms despite the large Coulomb barrier. How- 
ever, one can modify the Lagrangian by sending the weak 
boson masses to infinity to stabilize water. The resulting 
77 will be close to the experimentally measured value. It 
is not clear whether the same result can be obtained by 
simply discarding near ground states from the ensemble. 
Assuming this is correct, we investigate its consequence 
in our model by expanding <j> around the false vacuum 
above the T c of a first-order phase transition. As shown 
in Fig. 2, 77/s of the metastable state can be lower than 
the minimum 77/s of the true ground state (the effect will 
be more pronounced if c is bigger). The metastable state 
77/s has the tendency to keep decreasing until this state 
eventually disappears at higher T. It would be interest- 
ing to see whether this behavior persists in the strongly 
interacting theory. 
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